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Chapter  1 

INTRODUCTION 


Underwater  visual  imaging  is  often  restricted  to  very  short  range  by  suspended 
organic  material  and  particles  such  as  sand  and  silt  stirred  up  from  the  bottom.  In  this 
highly  turbid  environment  the  use  of  video  imaging  becomes  impossible  because  these 
suspensions  significantly  scatter  light.  Ultrasonic  waves  have  longer  wavelengths  and 
therefore  much  less  scatter.  The  necessity  for  divers  to  "see"  in  this  environment  for 
salvage  and  rescue  operations  and  military  applications  expresses  the  need  for  an  acoustic 
imaging  device. 

For  such  an  imaging  system  to  be  effective  several  factors  must  be  incorporated 
into  the  potential  design.  Real  time  operation  at  ranges  of  up  to  ten  meters,  good 
resolution,  and  a  large  field  of  view  are  all  desired.  The  equipment  must  be  lightweight 
out  of  water  and  as  close  to  neutrally  buoyant  as  possible.  The  design  should  be  as  simple 
as  possible  and  obviously  must  be  capable  of  operation  at  submergence  pressure. 

To  meet  the  demand  of  simplicity,  light  weight,  and  good  resolution,  an  acoustic 
thin  lens  and  a  zone  plate  will  be  developed  and  compared  for  use  as  a  beamformer.  A 
thin  lens  forms  beams  by  inserting  a  spatially  dependent  time  delay  as  the  wavefront 
passes  through.  Correctly  shaping  the  lens  will  insert  the  proper  time  delays  to  form  a 
clear  image.  The  lens  will  also  be  used  to  focus  the  transmit  beams  and  thus  eliminate 
most  of  the  expensive  and  complex  beamforming  electronics  of  the  imaging  system.  The 
Fresnel  zone  plate  uses  diffraction  principles  to  form  an  image.  This  paper  presents  a 
survey  of  the  design  parameters  of  thin  lenses  and  zone  plates. 
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Chapter  2 

THE  ACOUSTIC  THIN  LENS 


A  thin  lens  is  defined  as  one  that  has  thickness  very  small  compared  to  the  object 
or  image  distances.  There  are  several  parameters,  or  "degrees  of  freedom",  that  must  be 
determined  when  designing  a  lens  for  a  specific  application.  The  radii  of  curvature  of  the 
surfaces,  the  types  of  materials  to  use  and  their  indices  of  refraction,  the  diameter,  the 
thickness  desired,  and  the  focal  length  that  can  be  tolerated  are  all  parameters  that  may  be 
varied  to  obtain  the  desired  lens  performance. 

An  additional  requirement  for  this  application  is  size  and  weight.  It  is  always 
desirable  to  be  as  lightweight  and  as  small  as  possible.  The  complete  design  must  be  as 
close  to  neutrally  buoyant  as  possible.  This  puts  some  limitations  on  the  materials  that 
can  be  used. 


The  Basic  Thin  Lens  Formula 

The  formation  of  an  image  by  a  converging  lens  is  based  on  the  principle  that 
sound  will  obey  the  law  of  refraction,  Snell's  Law,  which  is  given  by: 

n,  sin0,  =  n2  sin02  ,  (1) 

where  n  is  the  index  of  refraction, 

0  is  the  angle  of  incidence  as  measured  from  the  lens  axis, 
and  the  subscripts  are  for  the  first  and  second  media. 

As  a  plane  wave  passes  through  a  converging  lens,  the  rays  are  bent  toward  the 
axis  of  the  lens.  At  some  distance,  /,  called  the  focal  length,  the  rays  will  converge  on  the 
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axis  at  a  focal  point.  Ideally,  if  the  plane  wave  is  placed  on  an  angle  0  to  the  axis,  the  rays 
will  focus  at  a  point  that  is  on  the  same  plane  (the  focal  plane)  as  the  original  focal  point 
perpendicular  to  the  lens'  axis. 


The  Lens  Makers'  Equation 

To  determine  the  radius  of  curvature  for  the  two  surfaces  required  to  produce  a 
desired  focal  length,  the  index  of  refraction  of  the  material  must  be  known.  These 
"degrees  of  freedom"  must  be  combined  into  a  form  that  can  be  used  to  design  and  create 
a  lens.  The  geometry  of  the  basic  thin  lens  problem  is  shown  in  Figure  1 .  The  refractive 
indices  of  the  three  materials  are  it\,  ni ,  and  n-i . 


Figure  1.  Basic  Thin  Lens  Geometry 
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A  ray  leaving  the  object  is  incident  on  the  first  surface  and  refracts  according 
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to  the  following  equation  for  spherical  surfaces  : 


n, 

s 


n2  +n, 
r 


(2) 


where  s  is  the  object  distance, 

s"  is  the  image  distance  to  the  center  of  the  lens, 

and  r  is  the  radius  of  the  spherical  surface. 

The  ray  then  strikes  the  second  surface  of  the  thin  lens  and  is  again  refracted  and 
directed  toward  a  new  image  point.  For  our  geometry,  the  thickness  of  the  lens  is 
negligible  compared  to  the  object  and  image  distances.  Therefore  combining  these 
equations,  we  obtain: 


(3) 


where  r\  is  the  radius  of  the  first  surface  and  rj  is  the  radius  of  the  second  surface. 

This  is  a  general  equation  that  the  lens  maker  can  use  for  a  thin  lens  with  different 
media  on  each  side.  This  equation  can  be  simplified  if  the  medium  on  both  sides  of  the 
lens  is  the  same: 
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When  the  surrounding  medium  is  water  (n,=l),  the  familiar  lens  makers'  formula 
is  obtained: 


_1_ 

/ 


=  («-!) 


1  1 


V/l  r2 


(5) 


where/is  the  focal  length. 


By  convention  ,  when  a  wave  is  traveling  to  the  right,  all  convex  surfaces 
encountered  are  considered  positive  radii  and  all  concave  surfaces  are  of  negative  radii. 


Ray  Tracing  Methods 

Graphical  ray  tracing  can  be  used  for  rapid  assessment  where  extreme  accuracy  is 
not  required.6  Basic  graphical  ray  tracing  is  done  by  drawing  concentric  circles,  shown  in 
figure  2a,  with  radii  proportional  to  the  indices  of  refraction  of  the  materials  used. 

Figure  2b  shows  an  example  of  graphical  ray  tracing  with  a  spherical  lens  surface. 


Figure  2.  Graphical  Ray  Tracing 
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In  figure  2b,  the  incident  ray,  a,  is  drawn  to  intercept  the  lens  surface  at  point  c.  In 
figure  2a,  a  line  is  drawn,  parallel  to  line  ‘a’  in  figure  2b,  from  the  center  of  the  concentric 
circles,  point  O,  to  the  circle  corresponding  to  the  medium  which  contains  the  incident 
ray,  labeled  A.  Next  a  line  parallel  to  the  radius  of  the  lens  surface,  line  c-d,  is  drawn 
through  point  A.  The  point  where  this  line  intersects  the  circle  of  the  second  medium  is 
labeled  B.  The  line  OB  in  figure  2a,  is  the  direction  of  the  refracted  ray,  b,  in  figure  2b. 
This  procedure  is  repeated  at  each  surface  until  the  ray  is  projected  onto  the  final  image 
plane.  This  is  a  good  first  check  of  the  system. 

When  greater  accuracy  is  required,  trigonometric  ray  tracing  is  used.  For 
spherical  surfaces,  the  path  of  a  meridional  ray  can  be  shown  quite  accurately.  This 
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method  is  used  in  several  computer  programs.  Kingslake  explains  this  procedure  in 
detail. 


Lens  Power 

The  power  of  a  thin  lens  is  the  sum  of  the  power  of  the  lens  surfaces.  The  power 
for  a  thin  lens  is: 


and  for  a  system  of  lenses. 


Power  =  —  =  (n  -  l)| 


Vri  r2j 


Power  =  T  — 

^  f 


(7) 


(6) 


7 


Lens  Aberrations 


Spherical  Aberration 

When  a  plane  wave  is  incident  on  a  thin  lens,  parallel  rays  can  be  traced  from  all 
points  on  that  lens.  It  can  be  shown  that  all  of  the  rays  do  not  focus  exactly  at  a  single 
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point.  This  image  problem  is  known  as  spherical  aberration  and  can  be  expressed  as: 

S  =  -h2/2R2[n^  f  ,  (8) 


where  ti\  and  «2  are  for  the  object  and  image  side  respectively, 
and  h  is  the  height  of  incidence. 


There  are  several  different  methods  to  keep  spherical  aberration  to  a  minimum. 
First,  the  shape  of  the  lens  can  be  altered  or  "bent".  Bending  the  lens  requires  that  the  lens 
power  and  therefore  the  focal  length  remains  constant.  From  equation  (6),  this  restriction 


requires  that 


1  1 


vri 


remains  constant. 


'2  J 


A  dimensionless  shape  parameter  is  then  defined  as: 


q  = 


h  +r\ 
r2  ~r\ 


(9) 


Equation  (9)  can  be  used  to  plot  shape  factor  vs.  axial  difference  in  focal  length. 
This  plot  is  typically  parabolic  in  shape  and  therefore  has  a  minimum  which  corresponds 
to  the  minimum  spherical  aberration.  Spherical  aberration  may  also  be  corrected  by  the 
use  of  multiple  surfaces  of  different  refraction  index. 
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Temperature  Dependence  of  Index  of  Refraction  (Thermal  Aberration) 

The  ultrasonic  index  of  refraction  is  the  ratio  of  the  speed  of  sound  in  water  to  the 
speed  of  sound  in  the  material.  Since  this  relationship  involves  the  speed  of  sound  it  is 
temperature  dependent.  For  optical  glasses,  this  temperature  coefficient  is  very  small,  but 
for  plastics  and  liquids,  it  can  be  quite  large.  This  can  cause  the  image  to  be  out  of  focus 
and  must  be  compensated  for.  The  thermal  aberration  has  been  treated  by  D.L.  Folds2  and 
is  summarized  below. 

Both  water  and  typical  plastic  lens  materials  have  relatively  large  values  of 
thermal  expansion  coefficient.  This  makes  the  lens  focal  length  highly  temperature 
dependent.  For  most  materials  the  sound  velocity  varies  linearly  in  the  temperature  range 
of  interest. 

The  sound  velocity  can  be  expressed  as: 


V  =Va+ IT  ,  (10) 

where  l  is  the  thermal  coefficient  of  sound  velocity  and  V0  is  the  velocity  at  0°  C. 

The  value  of  l  is  in  the  range  -1  to  -9  m/sec  °C  for  the  solids  and  about  4  m/sec  °C  for 
water.  The  index  of  refraction  for  the  material  can  then  be  expressed  by: 

n  =  n0  +bT  ,  (11) 


i4~lno) 


where 


b  =  n 


C, 


(12) 
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Combining  equations  (5)  and  (11),  the  variation  of  focal  length  with  temperature,  which 
is  our  point  of  concern,  is  then  expressed  by  : 


dI—*‘C.(n- 1) 


(13) 


The  axial  pressure  amplitude  response  for  circular  aperture  has  been  shown  to  be  the 
familiar  sine  function  relationship2,  specifically: 


Ar  =  20  log|sinc(7Tgx)|  ,  (14) 

where  g  =  (*  +  /)]  •  (15) 

In  this  equation  FN  is  the  lens  f/number  and  x  is  the  axial  distance  from  the  focal  point. 

Forcing  the  Strehl2  focusing  condition  that  Ar  remains  greater  than  0.5  dB  below 
the  maximum  response  level,  the  depth  of  focus  can  then  be  expressed  as: 


A F  = 


3  F„f 
D/A 


(16) 


Therefore  the  temperature  range  over  which  acceptable  focusing  will  exist  is: 


A  T  = 


A F 

f3F"l 

'cm-  1)' 

df/dT  ~ 

D/i) 

(4 -«„/)_ 

(17) 


This  equation  can  now  be  used  to  determine  an  acceptable  temperature  range  for  the  lens. 
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Thin  Lens  Performance 


The  measure  of  thin  lens  performance  is  its  directional  response  at  the  focal  plane 
and  its  depth  of  focus.  Axial  pressure  distribution  establishes  the  depth  of  focus  so  it  is 
common  to  plot  both  the  axial  pressure  distribution  and  the  focal  plane  pressure 
distribution.  For  the  current  application,  a  theoretical  thin  lens  was  designed  using 
polystyrene  and  the  following  prescription:  ri  =  -25  cm,  =  18.92  cm,  center  thickness 
of  1  cm  and  D  =  17.5  cm.  The  lens  equation,  equation  (5),  determines  a  focal  length  of  30 
cm  at  20°  C  for  this  lens.  Figure  3  shows  the  theoretical  directional  response  at  1  MHz  at 
the  focal  length.  The  beamwidth  is  near  the  0.8  degrees  (which  correlates  to  about  0.4  cm 
at  the  30  cm  focal  length)  expected  for  1  MHz  and  a  17.5  cm  diameter  lens. 


Figure  3.  Theoretical  Directional  Response  for  a  1  MHz  Polystyrene  Lens 


11 


The  other  key  parameter  in  lens  performance  is  depth  of  focus.  The  depth  of 
focus  is  predicted  to  be  about  3  cm  by  equation  (16).  Figure  4  shows  the  theoretical  axial 
pressure  distribution. 


Figure  4.  Pressure  Amplitude  (dB)  vs.  Axial  Distance  from  the  Lens  (cm) 
for  a  1  MHz  Polystyrene  Lens 
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Chapter  3 

THE  ACOUSTIC  ZONE  PLATE 


Introduction  to  Diffraction 

As  with  light,  sound  also  deviates  from  straight  line  propagation  when  passing  the 
edge  of  an  obstacle.  This  phenomenon  is  called  diffraction  and  can  be  explained  using 
wave  theory.  There  are  two  general  classifications  of  diffraction.  If  the  source  and 
observation  points  are  so  far  from  the  diffracting  aperture  that  lines  drawn  from  the 
source  or  observation  points  to  the  center  and  edge  of  the  aperture  do  not  differ  in  length 
by  more  than  a  small  fraction  of  a  wavelength,  then  it  can  be  treated  as  Fraunhofer 
diffraction.  If  this  is  not  the  case,  then  it  is  Fresnel  diffraction.11  While  there  is  no  exact 
wavelength  defined,  a  typical  cutoff  is  one  twentieth  of  a  wavelength. 

We  can  use  diffraction  principles  to  create  a  focusing  effect  using  a  Fresnel  zone 
plate,  shown  in  figure  5.  Each  zone  has  a  specific  radius  to  develop  this  focusing  effect 
using  diffraction.  A  Fresnel  zone  plate  was  designed  using  the  same  operating  frequency, 
1MHz,  as  the  thin  lens  designed  in  chapter  2.  It  contains  20  zones  and  has  a  diameter  of 
19.2  centimeters. 
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Figure  5.  A  1  MHz,  20  Zone,  Fresnel  Zone  Plate 


The  Zone  Plate  Governing  Equations 


Fraunhofer  Diffraction 

For  circular  aperture  with  normal  incidence,  the  Fraunhofer  diffraction  pattern  is 
axisymmetric.  If  the  variation  of  the  pattern  along  one  radius  of  the  observation  plane  can 
be  determined  then  the  rest  of  the  axial  patterns  are  known.  The  basic  geometry  is  shown 
in  figure  6. 
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ri  y  P 


Figure  6.  Basic  Geometry  for  Fraunhofer  Diffraction 

Since  we  have  determined  that  the  pattern  is  symmetrical,  we  need  only  determine 
the  pattern  along  the  a  axis.  It  can  be  shown11  that  the  amplitude  function  at  point  P, 
U0(P),  is  determined  by: 

U„  (P)  =  Ua  (oc0,P0  )Q(P)  ,  (18) 

where 

Q(P)  =  -^-T\\exp(-ikx  sin  6  2)dxdy  ,  (19) 

kD  j: 


and  the  double  integration  is  over  the  aperture,  a. 

Letting  x  be  constant  for  the  moment,  the  corresponding  limits  of  y  are: 


y  limit 


D 


—  X 


(20) 


and  the  integration  with  respect  to  y  yields: 

/<*= ■  (21) 

For  simplification  we  then  make  the  substitution: 

kD  . 

fi  =  — sin02 
2  K 


(22) 
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This  leads  to: 

4  \ 

Q(P)=^i 

/ 

K(d2  2  Y2  r 

2 - x  exp  - 

J4  J  L 

.2  K  1 
-i  —  iix  ax 

D  J 

If  we  let  the  integration  variable  be: 

il 

to 

(24) 

then 

r-Dx/ 

/2 

(25) 

dx  =  (D/2)dr  ,  (26) 

at  a  distance  of  x  =  ^  ,  the  integration  variable,  x  =  1 . 
Now, 

Q(P)  =  —  J(l  -  T2)^  exp[-i^ur}iT 


(23) 


(27) 


Dividing  the  integral  into  two  separate  integrals  from  -1  to  0,  and  from  0  to  1,  we  see  that 
the  value  of  ( 1  -  T2 )  is  the  same  for  both  and  the  sign  is  changed  in  the  exponential. 
Therefore  adding  the  exponentials  yields  2cos(/r/iT),  a  real  quantity.  Combining  and 
integrating  we  obtain  the  familiar  first  order  Bessel  function, 

Q(P)  =  n>0  .  (28) 

This  function  is  plotted  in  figure  7.  Since  this  pattern  is  symmetrical,  the  overall 
amplitude  pattern  is  a  set  of  concentric  rings  around  the  given  origin. 
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Figure  7.  Bessel  Function  of  First  Order 


The  amplitude  is  a  maximum  where  |i=0.  This  value,  |X,  is  related  to  radius.  If 
the  radius  is  r,  then  r  =  l  tan02  and  since  we  are  under  the  assumption  of  Fraunhofer 
diffraction,  we  have: 


tan  02  =  sin  02  , 
kD  .  . 

u  =  —  sin0,  , 
2k 


(29) 

(30) 

(31) 


Equation  (28)  was  first  derived  and  developed  by  G.B.  Airy  and  is  known  as  the 
Airy  disk.  Eighty  four  percent  of  the  total  energy  in  the  pattern  is  generated  within  the 
central  ring  where: 

1.22/A 


r  - 


D 


(32) 
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Fresnel  Diffraction 

With  the  need  for  a  short  focal  length  the  distance  from  the  diffracting  plane  to 
the  observation  point  must  be  necessarily  small.  The  assumption  that  the  rays  are  parallel 
from  the  aperture  to  the  focal  point  is  no  longer  valid.  Therefore  Fresnel  diffraction  must 
be  used.  Figure  8  shows  the  geometry  for  the  Fresnel  diffraction. 


T1  y  P 


Figure  8.  Geometry  for  Fresnel  Diffraction 

Since  the  diffracting  surface  is  circular  we  can  take  a  ring  shaped  element  2rcpdp, 
and  integrate  between  p=0  and  p=c.  This  leads  to  the  Kirchoff  integral  in  the  form: 

iA  c 

U0(P)  =  —  - — J  exp[il(r,  +  r2  )pnpdp 

*-ri0r20  0 


(33) 
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From  the  geometry, 

rt  =  rl  +  p 2  and  r2  =r2o+p2. 

Since  rlo  and  r2o  are  constants,  differentiating  and  dividing  equation  (33)  by  2  gives 

p  dp  =  r{drx  =  r2dr2 . 

Therefore, 

drx  =  fj/  jprfp  and  dr2  =  (j/  jp^p  , 

and  adding  we  get: 

pdp  =  -^—d(rx+r2)  .  (34) 

r,+r2 

Since  c  is  small,  we  can  make  the  approximation 

pdp  =  -^^—d(rx  +  r2)  .  (35) 

rio  +  r2o 


the  path  length  difference  being  expressed  in  half  wavelengths. 
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Changing  the  integration  variable, 


dl=AJ2dA , 


and 


since 


e'kl  _  e>k(r ,0+r2o)eirti 


kh=2n 


(38) 


(39) 


(40) 


The  first  exponential  is  constant  and  can  be  brought  outside  the  integral  and 
combined  to  form  U00(P),  the  amplitude  of  the  sound  with  no  diffracting  surface  present. 
Taking  A/2  outside  the  integral  we  get: 

A  (c) 

U0(P)  =  -ixU00(P)  Jexp(^A)JA  ,  (41) 

o 

We  have  A=0  at  1(0)  and  A(c)  is  when  p=c. 


Integrating  we  get: 

U0  (P)  =  Uao  (P){l-exp[/7rA(c)]}  ,  (42) 

and  the  corresponding  irradiance  is: 

I(P)  =  \uo(P)  |2  =  2/0 (P)[l - cos(^A(c))]  =  4/0 (P) sin2 A(c) j  .  (43) 


From  equation  (43),  it  is  readily  seen  that  I(P)  passes  through  a  series  of 
maximum  and  zeros  as  c  varies.  Therefore  the  pattern  is  a  series  of  dark  and  light  rings 
(fringes). 
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Fresnel  Zone  Plate 

Instead  of  taking  the  integral  over  the  whole  open  circle  in  equation  (41)  above, 
consider  a  zone  from  A=A,  to  A=A2: 

A2 

U0(P)  =  -iKU00(P)  Jexp(^A)JA  =  (/^[expO'flA,)  -exp(i7rA2)] 

Ai 

=  U00(P)[cos7tAl -costtA2  +/(sin^A, -sin^A2)]  .  (44) 

Now  let  A,  be  an  even  integer  and  Aj=A,  +1,  is  therefore  odd.  If  we  divide  the  complete 
circular  aperture  into  such  zones,  the  boundaries  of  which  have  A  equal  to  successive 
integers,  the  corresponding  radii  from  the  center  of  the  aperture  are  found  by  solving  the 
equation, 

A(c)  =  £tod  ,  (45) 

ho 


for  c,  with  A(c)  being  replaced  by  n  , 


n^or2o 

V  +  r20 


(46) 


If  we  cover  alternate  zones  with  an  sound  opaque  material,  the  result  is  a  Fresnel 
zone  plate.  To  be  able  to  compare  the  zone  plate  with  a  regular  lens,  if  we  make  r\0  the 
object  distance,  r2o  is  the  image  distance,  and  the  focal  length  of  the  zone  plate  is  given 
by 


/  rio  ho 


(47) 
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The  physical  meaning  of  focal  length  is  when  the  object  is  at  an  infinite  distance, 
the  image  is  formed  at  the  focal  length.  Using  this  in  equation  (46)  above  leads  to 

c\  =  nXf ,  or 


(48) 


This  is  known  as  a  geometric  zone  plate. 


An  interferometric  zone  plate  is  more  exact  but  more  difficult  in  practice.  As 
shown  above  ,  the  path  length  difference  to  the  focal  point  is  in  multiples  of  half 
wavelengths.  From  the  geometry  of  the  Fresnel  zone  plate,  shown  in  Figure  9,  and  the 

and  r2  =  F2  -  f 2 ,  therefore 


condition  above,  F  =  f  + 


nX 


(49) 


Typically,  f  »nX/  4  so  that  the  geometric  approximation  is  valid. 


Figure  9.  The  Fresnel  Zone  Plate  Geometry 
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Off  Axis  Aberration 


The  field  of  view  of  the  zone  plate  is  determined  by  the  off-axis  aberration.  The 
sound  path  difference  for  the  off-axis  case  is  shown  in  figure  10. 


Figure  10.  Off-axis  Aberration 


The  sound  path  difference  can  be  expressed  by: 


10 


rn  sin  a  +  f 


1>2 


1  +  ^tan  a  -  ry^j  -  /[l  +  tan 2  a]7^  (55) 


If  a  binomial  expansion  is  done  on  the  square  root  terms  and  the  sina  and  tana  are 
expressed  as  power  series,  the  sound  path  difference  becomes: 


rla 


8/ 3  2  f‘ 


3 rW_ 

4/ 


(56) 


where  the  higher  order  terms  have  been  dropped. 


From  equation  (48)  and  (49),  the  first  two  terms  are  equal  to  nX/2  so  the  third  and 
fourth  terms  are  the  off-axis  aberration.  As  before,  the  aberration  becomes  apparent  when 
the  sum  of  the  terms  becomes  A/4.  If  the  zone  plate  has  only  a  few  terms  (n  is  small),  the 
angle  can  be  made  large  and  only  the  fourth  term  will  be  significant.  Setting  this  term 
equal  to  A/4,  the  field  of  view  becomes 

a  =  {3nY'/2  .  (57) 


If  the  number  of  zones  is  large  the  third  term  dominates  and  the  field  of  view  becomes: 


a  = 


(58) 


Fresnel  Zone  Plate  Performance 

For  comparison,  a  theoretical  design  was  created  for  a  Fresnel  zone  plate  using  the 
same  criteria  as  the  thin  lens  of  Chapter  2.  With  a  focal  length  of  30  cm  and  a  frequency 
of  operation  of  1  MHz,  a  Fresnel  zone  plate  will  have  twenty  zones  with  the  radii 
determined  by  equation  (49).  Such  a  zone  plate  is  illustrated  in  Figure  5. 

As  with  the  thin  lens,  depth  of  focus  is  a  significant  concern.  Figure  1 1  shows  the 
theoretical  axial  pressure  distribution  for  the  designed  zone  plate.  The  peak  pressure  is  at 
the  focal  point  of  30  cm  as  expected.  The  first  side  lobes  are  down  at  least  15  dB  and  the 
depth  of  focus  at  the  3  dB  down  points  is  about  3  cm,  as  predicted.  Figure  12  is  the  focal 
plane  pressure  distribution  for  the  zone  plate.  The  beamwidth  at  the  3  dB  down  points  is 
approximately  0.2  cm  at  the  focal  point  and  the  first  side  lobes  are  down  about  18  dB. 
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Figure  11.  Axial  Pressure  Distribution  for  a  20  Zone,  19.2  cm  Diameter, 

1  MHz  Zone  Plate 


Figure  12.  Directional  Response  of  a  20  Zone,  19.2  cm  Diameter, 

1MHz  Zone  Plate 


Chapter  4 

SUMMARY  AND  COMPARISON 


26 


An  acoustic  lens  provides  an  inexpensive  way  to  focus  a  projected  beam  and  then 
focus  the  return  over  a  large  field  of  view.  The  images  are  not  degraded  by  suspended 
mud,  sand,  or  silt.  The  methods  of  optical  design  can  be  used  for  ultrasonic  lenses  and 
zone  plates  and  the  frequency  chosen  will  depend  on  the  range,  concentration  of 
suspended  particles,  and  the  resolution.  In  addition,  a  lens  with  low  spherical  aberration 
and  acceptable  temperature  insensitivity  can  be  designed  using  only  spherical  surfaces, 
which  are  easily  manufactured  to  close  tolerance. 

The  choices  of  materials  for  lenses  are  vast  and  some  inexpensive  and  easily 
manufactured  materials  have  been  used.1’2,3  Polystyrene  was  chosen  because  it  has  a 
good  impedance  match  with  seawater,  an  index  of  refraction  lower  than  water,  a  low 
insertion  loss,  and  good  physical  characteristics  for  a  lens.  Optical  lens  design  programs 
have  been  developed  and  can  be  used  to  refine  a  lens  design  for  ultrasonic  use.  Other  lens 
configurations  such  as  the  liquid  filled,  spherical  shell1  have  proven  to  be  excellent 
beamformers. 

The  zone  plate  is  a  little  more  difficult  from  a  manufacturing  standpoint  but  has 
some  of  the  qualities  necessary  for  a  good  imaging  system.  The  major  problem  will  be  to 
find  a  good  choice  for  a  material  that  is  “opaque”  to  sound.  Most  of  the  common 
anechoic  materials  have  some  transmission.  A  reasonable  design,  using  the  same 
parameters  as  the  lens  above,  was  developed  for  comparison  purposes. 

Figure  13  compares  the  theoretical  pressure  distributions  for  a  1  MHz  ideal 
polystyrene  lens  and  a  1  MHz  Fresnel  zone  plate.  The  thin  lens  is  a  slightly  better 
performer  than  the  zone  plate  as  it  has  a  slightly  narrower  beamwidth. 


Distance  along  the  Focal  Plane  (cm) 

Ideal  lens 
Zone  plate 


Figure  13.  Relative  Pressure  Response  of  a  1  MHz  Ideal  Thin  Lens 
and  Fresnel  Zone  Plate 
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Figure  14  compares  the  depth  of  focus  for  a  thin  lens  and  zone  plate.  Again  the 
thin  lens  is  a  slightly  better  performer  at  this  frequency.  The  zone  plate  may  have 
application  in  higher  frequency  applications,  such  as  medical  imaging,  where  the  zone 
plate  can  be  made  smaller  and  the  focal  lengths  are  much  shorter. 


0  10  20  30  40  50 

Distance  from  the  Lens  (cm) 

Zone  plate 
Thin  lens 


Figure  14.  Axial  Pressure  Distribution  for  a  1  Mhz  Ideal  Thin  Lens 

and  Fresnel  Zone  Plate 
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Appendix  A 
THIN  LENS  DESIGN 


This  Appendix  is  an  easy  to  use  "cookbook"  method  for  determination  of  approximate 
pressure  distributions  along  the  axis  and  in  the  focal  plane  of  a  thin  lens. 

The  pressure  distribution  along  the  axis  of  the  thin  lens  is  determined  as  follows: 

1  x  =- 19,- 18.9..  10  Define  the  distance  of  interest 

l(x)  :=f+-  x 

2.  g(X)  :=  - ^ -  Define  the  parameter  g(x) 

X---8-(fH-x) 

D 

3.  A(x)  :=  lsin(7t-g(x>'x)l  Define  the  formula  for  normalized  axial  pressure 

Jt-g(x)-x 


4.  P(  x)  :  =  20  log  |  ^  I  sirtfrr -g(  x)-x.)  [  j 1  £onvert  pressure  dB 


5.  Plot  the  axial  pressure  distribution 


20 


30 


40 


For  a  thin  lens,  the  pressure  distribution  in  the  focal  plane  is  determined  as  follows: 


1. 

x:=-2,- 1.9..  2 

2. 

c  :=  15000( 

cm/sec 

3. 

F:=l-106 

Hz 

4. 

X  :=— 

F 

5. 

D  :=  17.5 

cm 

f  =  30 

cm 

it  ' 

2-J1 - Dx 


4. 


P(x)  := 


U-f 


Jt 


•Dx 


X-f 


First  define  a  range  along  the  focal  plane 
Value  for  the  speed  of  sound 
Design  frequency 

Define  the  wavelength  you  are  using 

Define  the  lens  diameter  and  focal  length 

Equation  for  focal  plane  pressure 


5.  A(x)  :  =  201og(P(x)) 


Convert  to  decibels 


6.  Plot 


X 


Distance  along  focal  plane 


Appendix  B 

ZONE  PLATE  DESIGN 
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This  Appendix  is  an  easy  to  use  "cookbook"  method  for  determination  of  approximate 
pressure  distributions  along  the  axis  and  in  the  focal  plane  of  a  Fresnel  zone  plate. 

The  pressure  distribution  along  the  axis  of  the  zone  plate  is  determined  as  follows: 


1. 

2. 

to 

3. 


z=0,.1..40 

ua(z) 


:  =  ^/r(n 


s2  2 

)  +Z 


Define  a  desired  range 

Define  a  range  of  values  of  distance  from  the  radii 
the  axis 


ub(z) 


:  =  V  t(n 


.2  2 
)  tz 


Pj(z)  :  =  201og 


J^(j'sin(k-ub(z))  -  j*sin(kua(z))  -  cos(k-ub(z))  +-  cos(k-ua(z))) 
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Plot 


z 


For  the  zone  plate  the  pressure  distribution  is  determined  as  follows: 


1.  x=-2,- 1.9.. 2 


First  define  a  range  along  the  focal  plane 


2.  n  :=  1,3..  19  m(n):=n+l  Define  the  number  of  zones  (even  and  odd) 


3.  r(n)  := 


n-X-f+ 


n2-\2 


t(n)  := 


Aj 


m(n)-X-f-i- 


m(n)2-\2 


Determine  the  radius  of  each  zone 


4.  j:  =  V-l 


k  :  =  2-— 
X 


Define  wavenumber 


5.  R(n)  :  =  a/ f2  +  (r(n)  -  x) 


Define  a  distance  vector 


6.  P(x)  := 


2-7t 


rt(n) 


r(n) 


exp(jkR(n)) 

jk-R(n) 


•rdrdij)  Pressure  equation 


7.  A(x)  :=  20 log 


£exp(j-k-R(n))-(t(n)2-  r(n)2)- - * 

ritk-R 


(j(k-R(n))) 


.112637 


SPL  conversion 


8.  Plot 


Distance  along  the  Focal  Plane  (cm) 


